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Abstract 

Dissipative effects in electromagnetism on macroscopic scales are examined by coarse graining 
the microscopic Maxwell equations with respect to time. We illustrate a procedure to derive 
the dissipative effects on the macroscopic scale by using a Green-Kubo type expression in terms 
of the microscopic fluctuations and the correlations between them. The resulting macroscopic 
Maxwell equations are formulated within the General Equation for the Non-Equilibrium Reversible- 
Irreversible Coupling (GENERIC) framework, accounting also for inhomogeneous temperature. 
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1 



I. INTRODUCTION 



The dynamics of matter can be described on different levels of detail. For example, 

microscopic descriptions resolve short length scales and fast processes. In principle, one 

can follow the microscopic dynamics over very long times to arrive at the behavior on the 

macroscopically relevant time scales. However, doing so is often not desirable, and one 

can learn significantly more by using a coarse-grained description. For example, the flow 

of water in a complex geometry is usually studied not in terms of the dynamics of water 

molecules, but rather in terms of the continuum equations of hydrodynamics. This is so 

because key features of interest emerge only on the length and time scales much larger than 

the molecular ones. The intrinsic time scale for the relaxation of inhomogeneities in the 

macroscopic velocity field is closely related to the shear viscosity 77. This material specific 

property can either be determined experimentally on purely macroscopic grounds. Another, 

conceptually very interesting, route goes via the so-called fluctuation-dissipation theorem []]- 

4|, which establishes a relation between micro- and macroscopic properties. In particular, 

the Green-Kubo formula relates the shear viscosity to the correlations between fluctuations 

of the microscopic shear stress, a xy , by way of 
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V = j-f I ( or ^(^)° r ^(0)) eq rft, (1) 







with V the sample volume, T the absolute temperature, and (. . .) eq the equilibrium aver- 
age. The integration range [0, 00] stands symbolically for an integration range much longer 
than any microscopic relaxation time. However, it is implicitly assumed that this range is 
shorter than macroscopic time scales, which implies that micro- and macroscopic timescales 



y used in molecular dynamics sim- 
3j, where transport coefficients are 



are clearly separated. Relations of this type are frequent 
ulations to determine macroscopic transport properties 
expressed in terms of two-time correlations of current densities. 

The motivation for the work presented here related to the Maxwell equations rests on 
the following idea. The Green-Kubo relation ([T]), and generalizations thereof, not only 
tells us how to relate microscopic fluctuations to macroscopic transport coefficients. More 
fundamentally, it indicates that microscopic fluctuations give rise to dissipative processes 
in more coarse-grained descriptions. In view of electromagnetism, this serves as a strong 



motivation to believe that, e.g., fluctuations in the polarization and magnetization by way 
of particle vibrations and fast spin dynamics must be observable on macroscopic scales as 
dissipative processes. Although the Maxwell equations belong to the most studied differential 
equations in physics, they do not account for the dissipative phenomena beyond Ohmic 
resistance and the frequency dependent imaginary parts of the permittivity and permeability 
in the linear response regime. In view of the above discussion, we anticipate that the electric 
(Ohmic) resistance is related to the fluctuations of the electric current of the unbound 
charges. However, in addition also the bound particles fluctuate, giving rise to fluctuations 
in the polarization and the magnetization. How those fluctuations give rise to dissipative 
effects on the macroscopic scale is precisely what we want to elaborate on in this manuscript. 

Thermodynamics comes into play in the Green-Kubo relation ([I]) by way of the absolute 
temperature T. In addition, the occurrence of a "thermal" variable is also of fundamental 
importance in a purely macroscopic model with dissipative effects, since the latter lead 
to an entropy increase, and in turn to a change in the thermal state. In contrast, the 
Maxwell equations are usually taken as a set of isolated equations, decoupled from the other 
macroscopic variables, e.g., temperature and density. This means that the evolution of 
the electromagnetic field is generally obtained in the approximation that these macroscopic 
variables have set values or, at best, are given as functions of time. In order to obtain a 
rigorous and reliable theory, one needs to include these macroscopic quantities in the list 
of variables and to consider their dynamics in conjunction with the time evolution of the 
electromagnetic field. Only then one can capture the interplay between the thermodynamic 
behavior and the electromagnetic field, and dissipative effects in electromagnetism can be 
addressed in a consistent manner. 

To perform a thorough theoretical analysis of the interplay between the thermodynamic 
properties and the electromagnetic field is a demanding task, especially if one wants to 
account for dissipation effects beyond a linear response theory. Liu and coworkers have 
incorporated dissipative effects into the Maxwell equations 5Ml l| . The resulting theory 
describes the dynamics of macroscopic systems that are exposed to electromagnetic fields 
and contain electric charges and currents. In particular, the driving forces for the dissipative 
effects on the macroscopic scale have been established. In the dynamic equations, the new 
dissipative terms occur as an addition to the commonly used macroscopic Maxwell equations 
when the system is out of equilibrium. The validity of this theory was at first confined to 
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the low- frequency regime, but was later generalized to higher frequencies 11 1. 

In this paper the goal is to use the reversible Maxwell equations to derive macroscopic 
Maxwell equations for dissipative electromagnetism in a medium, valid also outside of the 
linear response regime. As the guideline to complete this task, the General Equation for the 
Non-Equilibrium Reversible- Irreversible Coupling (GENERIC) formalism of non-equilibrium 
thermodynamics 4 , [12I, [l3 ] is used for two reasons. First, thermodynamic consistency of the 
macroscopic model is ensured. And second, the formalism is equipped with a scheme for 
temporal coarse-graining, the most important part of which is concerned with the emergence 
of dissipative effects upon coarse-graining, i.e., a generalization of the Green-Kubo relation 
([T]). This scheme will be illustrated below specifically on the Maxwell equations in order to 
obtain the macroscopic equations that possess the GENERIC structure. These equations 
are then compared to the dissipative Maxwell equations suggested by Liu j^]. 

The paper is organized as follows. In Sec. Ull the microscopic and common macroscopic 
Maxwell equations are briefly presented. The usual coarse-graining procedure and its short- 
comings are recapitulated. In Sec. Ill II a GENERIC formulation of the Maxwell equations 
together with the equation for the energy density is presented. We illustrate the procedure 
to obtain, by temporal coarse-graining using Green-Kubo-type expressions, the dissipative 
effects in the time evolution equations for the macroscopic variables. 



II. MICRO- AND MACROSCOPIC MAXWELL EQUATIONS 

The classical derivation of the macroscopic Maxwell equations for the electromagnetic 
fields in a medium from the microscopic ones, 

b = -Vxe, V ■ b = , (2a) 

e = — V x b j , V ■ e = — p , (2b) 

is often performed through spatial averaging [l^] . Throughout the manuscript we denote the 
partial time derivative by A = J^A. The quantities eo and /io stand for the permittivity and 
the permeability, respectively, in vacuum. Starting from the microscopic Maxwell equations 
containing only fields e and b, one separates the free from the bound parts of the charge p 
and of the current density j. Upon spatial averaging, the expressions for the bound parts of p 
and j are simplified by introducing the polarization P and the magnetization M. With E and 
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B the spatially averaged electric field and magnetic induction, respectively, the macroscopic 
Maxwell equations can then be formulated conveniently as 

B = -V x E , V ■ B = , (3a) 

D = V x H - J e , V ■ D = p e , (3b) 

where one has introduced the electric displacement D = e E + P, assuming that quadrupol 
moments are negligible, and the magnetic field H = B/p — M. Doing so, the source terms 
of the macroscopic Maxwell equations are given by the averages of the free charge density 
and free current density, p e and J e , respectively. In addition, the evolution equation for the 
energy density of the electromagnetic field can be written in the form 

e em — — V • (E x H) — J e • E, (4) 

where E x H is the Poynting vector. 

The structure of the equations ([3]) imposes the interpretation that the field variables are 
D and B, the temporal equations are their time evolution equations, while the other two 
represent the constraints that must be satisfied at all times. The fields E and H then need 
to be expressed as some functions of D and B in order to render equations (J3J) closed. Such 
closing relations are known as constitutive equations, which contain the information about 
the medium, and are usually considered as given in macroscopic electrodynamics. The usual 
choice is to use linear relations between (E, B) and (D, H), valid if the electromagnetic field 
is sufficiently weak, in order to obtain the properties of the permittivity e and permeability 
fi that describe the response of a medium on the electromagnetic field. Then the real 
parts of the functions e(oj) and fi(oo) express the oscillatory motion of the bound charges 
and the imaginary parts describe dissipation. However, in many cases, like the passage of 
the laser beam through a substance, or for a medium in a strong magnetic field, the linear 
response theory is not valid anymore and one needs to work in the regime of nonlinear 
electrodynamics, with the consequence of loosing all the simple relations. 

A common constitutive relation concerns the current density. It is known that Ohmic 
conductors satisfy the condition E eq = at equilibrium, and that in non-equilibrium sit- 
uations the force to relax D to equilibrium is proportional to E. Such behavior is often 
expressed by J e = crE with conductivity a. This interrelation suggests to draw the following 
analogy. In stationary, equilibrium situations, with vanishing electric current J e , one obtains 
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from ([3]) the conditions 

V x H eq = , V x E cq = . (5) 

We therefore assume that in non-equilibrium situations on the macroscopic scale there are 
forces to drive B and D towards equilibrium with those forces being closely related to V x H 
and V x E with associated transport coefficients. Conditions (jSJ) have been obtained also by 
Liu by minimizing the total energy of the system under the constraints of the non-temporal 
Maxwell equations V ■ D = p e and V • B = [5(. 

The commonly used macroscopic Maxwell equations d3J) can be obtained not only by spa- 
tial averaging 
averaging 



ing 


see 


15, 


16| 



see e.g. |14|). which is also called "truncation" in |l5|, but also by ensemble 
161 ] . In both cases, only single time properties are considered, while tem- 
poral correlations and two-time ensemble averages are neglected. In the remainder of this 
manuscript, we will thus refer to these two techniques as "spatial" averaging, in contrast to 
temporal coarse-graining which will be discussed in detail below. Neglecting coarse-graining 
with respect to time to arrive at the spatially averaged Maxwell equations ([3]) may be a valu- 
able ansatz under many circumstances []jj]. However, it is also known that coarse-graining 
with respect to space and time can lead to a better understanding, and some aspects of a 

nn 

model emerge clearly only on longer time scales [4], 117| . Such an example is hydrodynamics, 
where the shear viscosity given by ([I]) and other transport coefficients in the Navier-Stokes 
equations arise at the continuum level of description, but are absent in a molecular model. 
Therefore, we will examine specifically the effect of coarse-graining in time by using Green- 
Kubo-type expressions to learn more about dissipative effects in electromagnetic systems 
with nonlinear constitutive relations. In order to obtain a closed set of the macroscopic 
time evolution equations for electromagnetism in a medium, a scheme for temporal coarse- 
graining of the "spatially averaged" Maxwell equations is presented following the GENERIC 



formalism 
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131 ] . The full coarse-graining procedure together with the expressions ob- 



tained for the additional dissipative terms is presented in the following sections. 



III. GENERIC FORMULATION 



A. Formalism 



The main points of the GENERIC framework of nonequilibrium thermodynamics 



6 



% ll2l, Ll3|] can be summarized briefly in the following way. In analogy to equilibrium thermo- 
dynamics, major importance comes to, first, choosing a complete set of variables x, which 
describes the situation of interest to the desired detail. The reversible contributions to the 
time evolution for this set of variables, x| rev , is formulated in close reference to classical 
Hamiltonian mechanics. In particular, it is related to the energy gradient by way of a Pois- 
son operator L, i.e., x| rev = L • (8E/8~x). The Poisson bracket associated to L, given by 
{A, B} = (<5v4/5x, L ■ 5B/5x) with appropriate scalar product (. . .), must be antisymmetric 
and satisfy the Jacobi identity, which are both abstract features that capture the nature of 
reversibility. In order to formulate the irreversible part of the dynamics one is motivated 
by the reversible part. For the reversible dynamics, the energy plays a distinct role. First, 
it is a conserved quantity for closed systems and, second, it drives the reversible dynamics. 
In parallel, for irreversible dynamics entropy is a fundamentally important quantity, which 
must not decrease for closed systems. In analogy to the reversible dynamics, it is assumed 
in the GENERIC framework that the irreversible contributions to the time evolution of x, 
x| irrev , is driven by the entropy gradient, i.e., that it is of the form x|i rre v = M- (5S/5x), with 
M a generalized friction matrix. The latter is required to be (Onsager-Casimir) symmetric, 
and contains transport coefficients and relaxation times associated to the corresponding dis- 
sipative effects. The condition that the friction matrix is positive semi-definite ensures that 
S > is fulfilled. 

In summary, the time evolution of x can be expressed in terms of four building blocks E, 
S, L and M as 

X = L(X ).^ + M (x).^>. (6) 

The two different contributions to the time evolution, reversible and irreversible, are not 
independent. Rather, they are interrelated by the two degeneracy requirements 

The first condition expresses the reversible nature of the L contribution to the dynamics, 
demonstrating the fact that the reversible dynamics captured in L does not affect the entropy 
functional. The second one expresses the conservation of the total energy of an isolated 
system by the irreversible contribution to the system dynamics captured in M. 

A particular feature of the GENERIC is that it is applicable on different levels of descrip- 
tion, including reversible Hamiltonian mechanics and dissipative macroscopic field theories. 



The corresponding four building blocks obviously differ between the different levels. How- 
ever, there are abstract procedures to relate them. Most importantly, the friction matrix 
on a coarse-grained level can be expressed in terms of the more microscopic dynamics ^]. 
This procedure helps to understand how dissipative effects can arise in coarse-grained de- 
scriptions. It is particularly this issue that we want to address in relation to the Maxwell 
equations of electrodynamics. 

B. Choice of variables, energy and entropy functionals 

First, we identify a natural set of variables for a dynamic description of dissipative elec- 
tromagnetism. In order to describe a body in an electromagnetic field, we make the local 
equilibrium assumption, so that temperature, internal energy and entropy densities are de- 
fined at each point of the nonequilibrium system. Although one is free to choose any of 
these three variables for describing the "thermal" state of the body, it will become clear 
below that it is particularly useful to choose the energy density e, since it is the density of 
a variable which is conserved for a closed system. The energy density e stands for the total 
energy density inside the system volume, consisting of both the internal energy of matter 
and the energy of the electromagnetic field. As far as the variables describing the electro- 
magnetic state are concerned, we keep in mind that our primary interest is in the effect 
of coarse-graining the Maxwell equations with respect to time. Hence, we use the already 
spatially averaged fields D and B as used in the temporal Maxwell equations in d3J). Since 
our focus here is on the temporal coarse-graining and on constructing the dissipative effects 
in a macroscopic formulation of electromagnetism, we consider only bodies at rest to avoid 
additional complexity, i.e., we do not include the velocity field in the set of variables. This 
restriction, in turn, also means that the mass density is a constant, as one can easily infer 
from the mass balance equation. In other words, we restrict our attention to materials for 
which the volume changes are negligible compared to the other effects of interest, namely 
materials with vanishingly small isothermal compressibility and thermal expansion coeffi- 
cient. The full set of variables is thus given by x = (e, D, B). The total energy and entropy 
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of the system, expressed in terms of these variables, are given by the functionals 



E 



S 




e d 3 r , 

s(e,D,B) d 3 



r 



(8b) 



(8a) 



The functional derivatives of E and S will be needed further below for the calculation of 
time evolution of variables x by using the general evolution equation (E]). In order to connect 
the expressions for these functional derivatives to known quantities, one refers to the Gibbs 
thermodynamic relation. The latter gives the change in the energy density e of a thermally 
isolated medium at rest in the form of the total differential 



under the assumption of local equilibrium, where the fields E and H, as well as the tem- 
perature T, are given as the partial derivatives of the energy density e with respect to the 
appropriate variables. The first term on the right side in is identical to the ordinary 
thermodynamic relation for the internal energy density in the absence of the electromagnetic 
field and for a constant mass density. The last two terms represent the change of the energy 
density due to changes in the electromagnetic fields. The form of these two terms is well 

nn 

known from the standard textbooks on electrodynamics [14], [18J, where the total electro- 
magnetic energy increment is given as a volume integral over the increments 5D and 5B 
multiplied by E and H, respectively. Relation ([9]) is also in agreement with the total energy 
density used in jjL9i] for a system at rest. We point out that the precise form of the energy 
density e of the system is not required below for the derivation of the macroscopic Maxwell 
equations in the GENERIC form. 

The functional derivatives 5E/5x. and SS/Sx are now obtained by reordering Eq.(j9]) in 
terms of the entropy total differential, and calculating the appropriate partial derivatives 
therefrom. One finds the following expressions needed for the GENERIC time evolution 



de = Tds + E • dD + H • d~B 



(9) 
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equation 



5E 
fa 



5S 
fa 



M 



V / 



(10a) 



1 

f 



I i \ 

-E 

V" H / 



(10b) 



C. Reversible dynamics 

The reversible contribution to the time evolution of the variables x = (e, D, B) is related 
to the energy gradient 5E/5x by way of the Poisson operator L. At this stage, we point out 
that the effects arising through coarse-graining in time by using the Green-Kubo formula 
are completely comprised in the friction matrix M, i.e., in the irreversible part of the time 
evolution. Correspondingly, temporal coarse-graining does not affect the reversible part of 
the time evolution of the variables x. Therefore, the Poisson operator 



V 



^ee L £ j} L £ b 

Lr> £ Ldd ^db 
L&e Lbd ^bb 



\ 



(11) 



can be build up from the already established Maxwell equations (J3J) and the adequate evo- 
lution equation for the total energy density, given only by the first term in Eq.flU) for a 
medium at rest. 

We approach the construction of L in the following manner, having in mind that the 
reversible contributions are of the form L • (5E/5x). Since only the first element of the 
vector 8E/5x. is non-zero, Eq. fllOap . terms L De and L Bs must reproduce the time evolution 
equations ([3]). The antisymmetry of L is then used to specify the elements L £ d and L £ b- 
With this, the e-component of the degeneracy requirement L ■ (5S/5x.) = can be satisfied 
by choosing L ee accordingly. Notice that this term is then also in consistence with the time 
evolution for e, and it is antisymmetric as well. The other two components of the degeneracy 
condition can be fulfilled by choosing Ldb = Vx, Led = -Vx, and Ldd = ^bb = 0. In 
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summary, one obtains the Poisson operator 



V 



E- V x H-H- V x E 
V x H 
-V x E 



-H-Vx E-Vx 
Vx 
-Vx 



(12) 



with all derivative operators acting on everything to the right, i.e., also on functions mul- 
tiplied to the right side of the operator L. At this point we mention that, in general, the 
symbol "•" in (E]) implies not only summation over discrete indices. If field variables are 
involved the operators L and M are written in terms of two space arguments (r, r'), and an 
integration over r' must be performed when multiplied with a function of r' from the right. 
However, in the case of the field equations being local, one can express L and M in terms of 
a single variable r only [4], and no integration is implied when these operators are multiplied 
from the right. Such single variable notation is used for the Poisson operator f[T2"j) . However, 
in case of the friction matrix M, discussed in the next section, it will be beneficial to keep 
the general form in terms of r and r'. 

The reversible time evolution of the variables x as obtained from L ■ 5E/x takes the form 



i rev = -V • (E x H) , (13a) 
D rev = V x H , (13b) 
B rev = -V x E , (13c) 

in accord with ([3]). 

We note here that the lower-right 2x2 block matrix in L, Eq.f ll2p . describing the 
electromagnetic part of the system, has been used previously by Marsden and Ratiu in 



2CH ] . Therein, it was given in terms of a Poisson bracket at which we can arrive by using 
{F, G} = j(5F/5x) ■ L • (5G/5x.)d 3 r. In our work, their Poisson bracket is embedded in a 
thermodynamic formulation with a thermal variable e. 

The degeneracy requirement and the antisymmetry of the Poisson operator are satisfied 
by construction. For the Jacobi identity, it is useful to observe that in our system the 
time evolution of the entropy density s (e, D, B) has no reversible contributions. Since the 
transformation of variables does not affect the Jacobi identity, we choose to prove the Jacobi 
identity in the set of variables x' = (s, D, B), in which the Poisson operator L' takes a much 
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simpler form. Namely, only two elements, L'db and L'bd, are non-zero, and L' coincides 
with the Poisson bracket in 20|, for which the Jacobi identity is fulfilled. With this, the 
Poisson operator L is completely constructed. 



D. Friction matrix obtained from microscopic fluctuations 

Coarse graining with respect to time from one level of description (LI) to another one (L2) 
in principle leads to two kinds of irreversible effects on the coarser level. First, dissipative 
effects on level LI reappear on level L2. Second, certain effects that are slower than the 
time resolution on level LI and faster than the time resolution of L2 will be expressed as 
rapid fluctuations on level L2, and hence emerge as irreversible effects on that coarser level 
of description. It is this latter kind of irreversible effects after coarse-graining which we 
concentrate on. The coarse-graining procedure is captured in the following five steps and 
illustrated for the case of dissipative electromagnetism. 

First, we express the friction matrix M, which captures the dissipative contributions 
to the GENERIC evolution equation ([6]), in terms of a Green-Kubo-type expression. The 
matrix M can be calculated from the microscopic description by a generalization of the 
Green-Kubo relation for the viscosity (Q. If x f denotes the fluctuations on all components 
of x, we write jj], 



M(r,r') = / T rft(x f (r,t)x f ' T (r',0)) , (14) 

fc B Jo 

which can be expressed alternatively, in component notation, as 

[l + e(x i )e(x k )]M ik (r,r') = -^-(Axl(v)Ax[ (r')> , (15) 

with e(xj) = ±1 depending on the parity of Xj under time- reversal. Here, r is an intermediate 
time scale separating the slow degrees of freedom (on the coarser level L2) from the fast ones 
(on the finer level LI). The quantity AxJ(r) is given by AxJr) = f Q x~(r, t) dt, and the 
brackets (. . .) indicate the average over an ensemble of microscopic trajectories consistent 
with the slow macrostate x over the whole time interval, see pages 357-358 in [4] for more 
details. The superscript T in ffT4"]) denotes the usual matrix transpose. 

Second, for formulating the form of the fluctuations x f the type of the variables x is 
important. Specifically, choosing densities of (conserved) extensive variables is particularly 
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useful. Doing so is in close ana 
in epical ther.odyna.nics Q 



ogy to the common procedure in the theory of fluctuations 
2ll | . Furthermore, the Green-Kubo relation ([I]), and similar 
relations for the thermal conductivity and diffusion coefficient, indicate that current-current 
correlations are quantities of fundamental importance. For example, the stress tensor is the 
current of momentum. Current densities play a major role in the formulation of conservation 
laws in the form of field equations for density variables of extensive quantities, e.g., internal 
energy density, momentum density. Hence, from this perspective it is tempting to use 
density variables in the set x. How this can be exploited further is discussed in the third 
step below. In our example, the energy density e is hence a good variable, in contrast to 
the temperature T, or the (not conserved) entropy density s. Further dissipative processes 
we are considering here originate from the fluctuations in the particle positions and spins. 
Thus, fluctuations in D and B arise due to fluctuations in densities of extensive variables, 
namely in the polarization P and the magnetization M. 

Third, we make an ansatz for the structural form of the fluctuations. We assume that 
the evolution equations for the fluctuations are similar in structure to their macroscopic 
counterparts, in the spirit of Onsager's regression hypothesis. This hypothesis states that 
the fluctuations about the equilibrium state decay, on the average, according to the same 
laws that govern the decay of macroscopic deviations from equilibrium 2|-|4|. In view of 
the first component of x f we note that (j3J) represents the change in electromagnetic energy 
only, while the corresponding equation for the thermal energy, e t h, of the body will be of 
the form e t h = J e • E — V • J q with heat flux 3 q , such that the total energy, i.e. the integral 
of e — e cm + e t h, is conserved. Approximating the fluctuating contributions by white noise 
leads to expressions for Ax f of the form 

1 _v ■ (W( q ) + E x WW + W( E ) x H) ^ 
V x W< H > - W® 



Ax 1 



y -V x W( E ) j 



(16) 



where all W are (small) increments of Wiener processes, which are integrals of white noise 
over a time interval r. The superscripts to the Wiener processes indicate their physical 
origin. 

Fourth, by specifying the second moments of the Wiener processes, and all correlations 
between them, the construction of the M matrix is completed. We assume that all W are 
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decorrelated because of the difference in the underlying microscopic mechanisms, but that 
W q has a non-zero correlation only with the electric current fluctuations W-i, which will 
give rise to the thermoelectric effects. Hence we write the following relations 



(W^(r)W^(r')) = 2A; B r5(r-r')C (a/3) (r) 



(17) 



with a,/3 G {q,j,E,H}, where the non-zero i.e., C( qq \ C®, C< EE \ C< HH > and 

= C ( --' q ^ T , represent the correlation functions, with implicit r-dependence, e.g., through 
temperature. Note that we have assumed only local correlations between the Wiener pro- 
cesses, i.e., between the fluctuations. Hence, we are dealing with a local field formulation. 
However, non-local effects can be included if so desired. 

In the fifth and final step, the firction matrix M is calculated. Using the fluctuations 
(TiT)]) in conjunction with the correlations ffTTl) we calculate the Green-Kubo relations for the 
elements of the matrix M(r, r') by way of (fl5l) . The final result for the friction matrix is 

/ /i/f /i/f _ fl/f _ \ 



M(r,r') 



M ££ 


M eD 




M De 


M DD 





M Be 





M BB 



(18a) 



with 



M ££ 

M DD 
M bb 

M eD 
M De 
M eB 
M Be 



V ■ C (qq) (r) ■ (V'5(r - r')) + V ■ E(r) x C (HH) (r) ■ (V x E(r')5(r - r')) 

+V ■ H(r) x C (EE) (r) ■ (V x H(r')5(r - r')) , (18b) 

-V x C (HH) (r) • (V8(r - r')) x + 5(r - r')C (jj) (r) , (18c) 

-V x C (EE) (r) ■ (V5(r - r')) x , (18d) 

V • C (qj) (r)(5(r - r') + V • E(r) x C (HH) (r) • (V'<J(r - r')) x , (18e) 
C (qj) - T (r) • (V'5(r - r')) - V x C (HH) (r) • (V x E(r')5(r - r')) , (18f) 

V • H(r) x C (EE) (r) ■ (V'5(r - r')) x , (18g) 
-V x C (EE) (r) ■ (V' x H(r')5(r - r')) , (18h) 



where it is understood that contractions • and cross products x are executed in the order of 
occurrence from the right to the left. 

Given the expression (fl8"|) for M, it is straightforward to show that the degeneracy con- 
dition M • (5E/5yi) = is indeed satisfied. The origin of this result goes back to the ansatz 
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for the fluctuations Ax f in (|T6j) . Their specific form ensures that the fluctuations occur on 
a submanifold of constant energy, as can be shown readily by 

Ax f - T .f A = 0, (19) 
ox 

as required by the GENERIC. Equation f lT9|) is fulfilled already by writing the first compo- 
nent of x f , namely the fluctuations in the total energy density, as in (Tl6l) . which requires it 
to be the divergence of a vector field. It is then straightforward to show that ( 1T9|) is satisfied 
if boundary terms can be neglected. 

If dissipative matrix M(r, r') is applied to the entropy gradient 5S/5x(r'), including an 
integration over r', one gets the irreversible contributions to the time evolution equations. 
Using the definitions 



E* 
H* 



ic< EE ) ■ (V x H) , (20a) 

_I C ( HH ) • (V x E) , (20b) 

^ = COO ■ - + C^> T • vi , (20c) 

one obtains for the irreversible part of the time evolutions of D and B 

D irr = Vxff-j;, (21a) 

B irr = -V x E*, (21b) 

where all functions are evaluated at r. We recognize in D; rr the Ohmic conduction C^, 
while represents the thermoelectric coupling |l|. 



E. Final set of evolution equations 

The final electromagnetic time evolution equations can be conveniently written in the 
following way. Combining the reversible (fT3l) and irreversible (I2T]) contributions to the time 
evolution equations of the fields D and B, the resulting temporally coarse-grained Maxwell 
equations differ from ([3]) in that the fields E, H and J e in ([3]) have to be replaced by E + E*, 
H + H* and J e + J*, respectively. The additional fields E*, H* and J* give rise to the 
irreversible contributions to D and B. They are proportional to the quantities V x H, 
V x E, E and VT, as expected from conditions (J5]) and the discussion in Sec. [TTl and 
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represent the nonequilibrium forces that tend to restore equilibrium. We mention that the 
full macroscopic Maxwell equations obtained here can be used to derive expressions for the 
complex permittivity e(u) and permeability when examined in the linear response 

regime 0, Q] . 

The frequently considered balance equation for the energy density can be obtained from 
(I13al) and the result of the application of M(r, r') to the entropy gradient 5S/5x(r'), namely, 



-V • (E x H) - V • x H* + E* x H + C (qq) • + C (qj) • . 



(22) 



The first term on the right side, which is the divergence of the well known Poynting vec- 
tor, represents the reversible part of the time evolution equation, while the second term is 
the irreversible part, describing the heat conduction (through C^ qq ^) and the Peltier effect 
(through C^), as well as a modification of the electromagnetic energy flux, i.e., of the 
Poynting vector, due to the additional fields E* and H*. This modification of the energy 
flux was found also by Liu [7, 9]. 

For the time evolution of the entropy density s one can write, using the chain rule and 
rearranging the expression, 



v . ( Ic( qq ) • v- + -c^) • -) 

\T T T TJ 

/ vl \ / c (qq) c (qj) 



+ 



E 

T 



|(VxE) 
V i (V x H) ) 





o o c( HH ) 









\ 



/ 



\ oo c( EE ) / 



E 

r 



\ 



|(VxE) 
V i (V x H) J 



(23) 



The central matrix of the last term is positive semidefinite as can be seen after inserting 
the microscopic expressions for the correlation functions ( jTTl) . and therefore the entropy 
production rate is indeed non-negative. The obtained entropy production rate corresponds 
to the one proposed by Li 



I non-negati 
jfu in {^Jf]]. 



IV. DISCUSSION AND CONCLUSIONS 



Coarse-grained dissipative Maxwell equations of electromagnetism have been obtained 
and examined from the perspective of nonequilibrium thermodynamics by using the 
GENERIC formalism. We have illustrated a scheme to relate the dissipative effects on 
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the macroscopic scale to the fast microscopic fluctuations in the polarization and magneti- 
zation in the spirit of a generalized fluctuation-dissipation theorem. In particular, the study 
is applicable also to materials having nonlinear constitutive relations between (E, B) and 
(D, H). It emerged clearly that the difference between spatial and temporal coarse-graining 
is of fundamental importance. On the one hand, the usual procedure of averaging the mi- 
croscopic Maxwell equations only spatially (or by single-time ensembles averages) leads to 
terms equivalent to the reversible contributions in the GENERIC formulation, L ■ (5E/5x). 
On the other hand, coarse-graining the spatially averaged Maxwell equations further with 
respect to time leads to new irreversible effects captured in terms of M • (SS/Sx), as shown 
here. The friction matrix M has been related to two-time correlations of fast processes 
on the finer level of description. Using the Green-Kubo expression (|14|) . the correlations 
of microscopic fluctuations give rise to dissipative processes such as Ohmic currents, the 
thermoelectric effect, and to other irreversible contributions to the electric and magnetic 
fields. 

In the above investigations, the following points have to be highlighted. First, we pre- 
sented an illustration of the coarse-graining procedure in order to obtain the structure of 
fluctuations and, in turn, the structure of the dissipative effects in electromagnetism. Sec- 
ond, the identification of current densities occurring in x f was of fundamental importance, 
since correlations between them lead to the new dissipative processes on the macroscopic 
scale. To aid the relation between the Green-Kubo-type expression for M ( fT4l) and the 
correlations between current densities, it has proven useful to choose density variables of 
extensive quantities in the set x. Using Onsager's regression hypothesis, the evolution equa- 
tions of the fluctuations takes the form of conservation laws, including fluctuating current 
densities. We point out that the correlations f }TT|) are similar in spirit to the Green-Kubo 
relation ([1]), in the sense that they define transport coefficients. In contrast, the expression 
for M ( [14]) in conjunction with the entropy gradient leads to the full form of the dissipa- 
tive processes. The third point concerns the details behind the correlations ( TlTl) . In order 
to specify the dissipative processes beyond their structure, one needs to discuss in more 
detail the physics behind the correlation matrices C a/3 , with G {q,j,E, H}. That is, 
the transport coefficients appearing in the final expressions through C al3 can be obtained, 
e.g., via molecular dynamics simulations, or be derived based on explicit expressions for the 

—pie fluxes and the —o. between the. H. 
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The scheme of temporal coarse-graining for obtaining the macroscopic Maxwell equations 
differs from Liu's (5), [g], who started by making an ansatzioi the entropy production rate. 
However, both procedures lead to identical modifications of the Maxwell equations, and also 
the energy equation (I22p is in agreement with the corresponding relations in The 
obtained Maxwell relations and extensions of them are successfully applied to the problems 
of nematic liquid crystals (5), Q], colloidal magnetic and electric fluid, and ferrofluids (3- 
9[. In nematic liquid crystals, the dissipative couplings between the new thermodynamic 
forces V x E and V x H with the additional thermodynamic forces specific for this system 
are permitted and are equipped with transport coefficients, thus being valid also in case of 
the oscillatory instabilities, when usually employed static Maxwell equations no longer hold 
Liu's equations for ferrofluids also give rise to the dissipative forces that account for 



the special spin-up behavior of ferrofluids under a rotating external field |t - Q] . 

A cross relationship, similar to the one suggested by Liu (5)], occurred in the dissipative 
terms ( !2~T1) obtained by us: C^ HH ^ in Di rr accounts for the correlations of the fluctuations 
in magnetization, i.e. the magnetization relaxation time, and C( EE ) in Bj rr accounts for the 
correlations of the fluctuations in polarization, i.e. the polarization relaxation time. This is 
the consequence of the nature of Maxwell equations and of the mutual dependence of the 
electric and magnetic fields. It is naturally obtained by the form ( TT6|) for the increments in 
fluxes Ax f . 



Aspects, different from the ones discussed here, of dissipative effects in e 



ectromagnetism 



23| . Decomposi- 



and of linking electromagnetism with hydrodynamics can be found in 
tion of the electromagnetic fields, in particular polarization and magnetization, into the slow 
and fast parts has been considered by Mderhof and K r oh Q. ,n this way they extended 
the irreversible thermodynamics approach used by de Groot and Mazur jl| for the hydrody- 
namics of magnetic and dielectric fluids in interaction with the electromagnetic field. In a 



subsequent paper [231 ] . Felderhof studied semirelativistic hydrodynamic evolution equations 



for a medium with polarization and magnetization. 
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